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_ ■ We consider rectangular graph superlattices of sides l\,l2 with the wavefunc- 

dj . tion coupling at the junctions either of the 5 type, when they are continuous 

, and the sum of their derivatives is proportional to the common value at the 

| junction with a coupling constant a , or the 5' s type with the roles of func- 

tions and derivatives reversed; the latter corresponds to the situations where 
■ the junctions are realized by complicated geometric scatterers. We show that 

. the band spectra have a hidden fractal structure with respect to the ratio 

Q := t\l li . If the latter is an irrational badly approximable by rationals, 5 

ON 
O 



lattices have no gaps in the weak-coupling case. We show that there is a 
I/-) ' quantization for the asymptotic critical values of a at which new gap series 

open, and explain it in terms of number-theoretic properties of 6 . We also 
show how the irregularity is manifested in terms of Fermi-surface dependence 
on energy, and possible localization properties under influence of an external 
_L ! electric field. 

£ 

O 
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: I Introduction 

The advent of new technologies which made possible fabrication of semiconductor 
quantum wires and other tiny structures opened a new chapter in solid-state physics 
as well as ways to many potentially useful devices. At the same time, this develop- 
ment has an impact on the quantum theory itself, which is maybe not so spectacular 
but by no means less important. The point is that by investigating various "tay- 
lored" systems one is able to study — both theoretically and experimentally — 
interesting and sometimes unsuspected effects "hidden" in the basic equations of 
quantum mechanics. 

An example is represented by irregular spectral properties due to incommensu- 
rability of certain parameters of the system. Such behavior is known to occur, e.g., 
for two-dimensional lattice electrons in a constant magnetic field. 1-5 In this paper 
we are going to discuss similar effects which can be observed without presence of an 
external field in certain graph superlattices. 

Before we shall describe the model, a few words should be said about simplifi- 
cation it involves. A real quantum wire is a complicated many-body system; even 
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Figure 1: A rectangular lattice 



if we suppose that it has an ideal crystallic structure, we may describe it as an elec- 
tron duct at most in a certain range of energies where the profile of the conduction 
band is reasonably flat. Neglecting the lateral size of the wire, i.e., assuming that 
the propagating electrons remain in a single transverse mode, is another approxi- 
mation. It can be justified in thin wires, not only by practical experience, but also 
by rigorous arguments 6 showing that the intermode coupling involves a dynamical 
(p-space) tunneling, and therefore it diminishes exponentially with the decreasing 
wire thickness. 

Although the replacement of a quantum wire system by the corresponding graph 
structure leaves us with a much simpler model, other idealizations may be useful 
to draw lessons of it. A typical one concerns the global size of the system. As the 
number of individual cells in a superlattice grows, "collective effects" become more 
important; the question is when they become to prevail. The larger the superlat- 
tice, the more reasonable is to start from an infinitely extended structure, with the 
boundary effects considered perturbation. 

The last introductory remark concerns the question why various results based on 
number-theory properties of the parameters represent more than a nice mathemat- 
ics. With a certain exaggeration one can certainly claim that for a mathematician 
all rationals are the same, while in physics with its finite-resolution experiments it 
is meaningless to ask whether a measured quantity is irrational or not. Fortunately, 
differences between number types and between simple and complicated rationals are 
usually two sides of the same coin. 

While theoretically there might be an ultimate difference between the value 
\/2 and its close rational approximations, it is manifested over a long scale (in 
energy, time, etc.); at a shorter scale it is important that the value in question does 
not coincide with one of simple harmonies. The irrational model thus represents a 
proper description of the "disharmonic" situation. These somehow vague statements 
can be readily illustrated with the help of the system treated below. 

II Description of the model 

The configuration space of our model is a two-dimensional lattice graph whose 
elementary cell is a rectangle of sides £i, £2 {cf. Figure 1). If no external field is 
applied, the motion of electrons on graph links is free. Since the choice of the energy 
scale will play no role in the following, we choose atomic units, h 2 /2m* = 1 , for 
the sake of simplicity; hence if the wavefunction ip is supported in the interior of a 
single graph link, the Hamiltonian changes it into — ip" . 
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The nontrivial part of the problem concerns, of course, the behavior at the junc- 
tions. The wavefunctions must be coupled there in such a way that the probability 
flow is preserved; in mathematical terms this is equivalent to the claim that the 
Hamiltonian is a self-adjoint operator. This is known for long: Schrodinger opera- 
tors on graphs appeared in quantum mechanics for the first time in connection with 
the free-electron model of organic molecules, 7 and in recent years interest to them 
has been revived. 8 " 16 

The requirement of probability-flow conservation does not specify the coupling 
uniquely: it can be satisfied, e.g., if the wavefunctions are continuous at all vertices 
and their derivatives satisfy there the conditions 



where m is the vertex number, the derivatives are taken in the same direction 
(conventionally outward), the sum runs over all links entering this vertex, ip(x m ) is 
the common value of the functions ipj there, and the real numbers a m are coupling 
constants characterizing the junctions. For the sake of brevity, we shall refer to 
(Refdelta be) as to the 5 coupling; the name is motivated by the fact that in the 
simplest case of just two links this is nothing else that the 5 interaction on a line. 2 

A choice of the coupling should, of course, be obtained by deriving the graph 
model from a more realistic description, in which the configuration space consists of 
a system of coupled tubes. Though a heuristic argument showing that for an ideal 
star-like junction the conditions (Refdelta be) with a m = might the optimal 
choice was given more than four decades ago, 7 and it is natural to expect that 
nonzero coupling constants could correspond to local deformation of the junction 
region, impurities or influence of external fields — in short, imperfect contacts — no 
convicing answer is known up to now. 

Moreover, the coupling (Refdelta be) is not the only possible: for a junction of 
N links there is in general an iV 2 parameter family of self-adjoint operators which 
act as the free Hamiltonian outside the branching points. A method to construct 
such operators and some particular classes of them were discussed in detail in Ref.14. 
In distinction to (Refdelta be), the "additional" couplings have wavefunctions dis- 
continuous at the vertex, i.e., the limits for at least one pair of links differ mutually 
there. 

This feature is not automatically disqualifying. It was shown in Ref.10 that 
the so-called 5' interaction 17 represents a reasonable (if idealized) model for a com- 
plicated geometric scatterer, in which instead by a point contact two halflines are 
joined by numerous short "wires". Moreover, this result extends to junctions with 
any number of links; 13 the most natural counterpart to (Refdelta be) appears to 
be the so-called 5' s interaction which requires the wavefunction derivatives to be 
continuous, tjj[{x m ) = ■ ■■ = ifj' N (x m ) =: ip'{x m ) , and 
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for some (3 m . The "coupling constants" (3 m here measure, roughly speaking, the 
total length of the wires which constitute the geometric scatterer; for a more detailed 
discussion we refer to Ref.13. 

In what follows we shall be concerned with rectangular lattices in which the 
coupling at each junction is the same and belongs to one the above described types; 
for the sake of brevity we shall refer to them as to the 5 and S' s lattices, respectively. 
In a sense, such lattices represent a generalization of the classical Kronig-Penney 
model and its 5' modification 2 to higher dimensions. 18 



Ill General properties of the spectra 

Before proceeding further let us recall some results about the spectra of the consid- 
ered lattice Hamiltonians derived in Ref.13. 



III.l 5 lattices 

By assumption, a 5 lattice is a periodic system in both directions. Performing the 
Bloch analysis, we arrive at the band condition 19 

cos -&i£i — cos k£i cos $2-^2 — cos k£ 2 a , , 
_ £ _i_ _ £ _ — n (i\ 

sinA;^ sinH 2 2k K } 

Although an analytic solution can be written in the trivial case only, the condition 
(Re fdelta band condition) nevertheless allows to draw many conclusions about the 
spectrum. Let us rewrite it in the form 

a JU Vj — cos k£j 
2k ~ p{ sin k£j ' 

if the quasimomentum components -&j , j = 1, 2 , run though the Brillouin zone, the 
ranges of the parameters Vj := cos-d j£j cover the interval [—1, 1] . It is easy to see 
that for a fixed k , the maximum of the right side equals 



^ / * \ x — > ( k£j 7i 
F+(k) := J2 tan 
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where the the square bracket denotes conventionally the integer part, and the mini- 
mum, F-(k) , is given by a similar formula with tan replaced by — cot . It is clear 
from here that the gaps of the 5-lattice spectrum with a coupling constant a on 
the positive part of the energy axis are determined by the condition 

± Y k > ±F±{k) (4) 

for rka > , respectively. The negative part of the spectrum is obtained analo- 
gously by comparing y~ with the extremum values of the function iF±(in) . Simple 
consequences of the condition (Re fdelta gap) are the following: 
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(al) The spectrum has a band structure; it coincides with the positive halnine 
[0, oo) if and only if a — . 



from the second one, equals k 2 n or k 2 m . 
(a3) The lowest band threshold is positive for a > and negative if a > ; in 



(a4) The positive bands shrink with increasing \a\ . 

(a5) Each gap is contained in the intersection of a pair of gaps of the Kronig-Penney 
model with the coupling constant a and spacings t\ and £ 2 , respectively. 

(a6) All gaps above the threshold are finite. If there are infinitely many of them, 
their widths are asymptotically bounded by 2\a\(£ 1 +£ 2 )~ 1 + 0(r~ l ) , where r 
is the gap number. 

Most of these results have a natural meaning. In particular, (a5) shows that trans- 
port properties of the lattice are better than a combination of its one-dimensional 
projections. Notice that the Kronig-Penney spectral condition 2 can be cast into the 
form (Refdelta band condition) with a single trigonometric expression on the Ihs. If 
an energy value is contained in a Kronig-Penney band in one of the directions, it is 
trivially also in a band of the lattice Hamiltonian, the other factor being annulated 
by choosing 9j — k . The converse is not true, of course: the condition (Refdelta 
band condition) may be satisfied even none of the factors can be annulated sepa- 
rately. The directions in which the electron is able to "dribble" through the lattice 
will be seen in Section VI below. 

Less trivial is the irregular dependence of the spectrum on the rectangle-side 
ratio 9 := coming from the existence of competing periods in F±(k) . It 

appears that it is not only rationality or irrationality of 9 which matters, but also 
the type of irrationality plays a role. In this respect, the situation is similar to the 
almost Mathieu equation mentioned in the introduction. 

Let us recall some elementary facts from the number theory. 20 An irrational 
number 9 is badly approximable by rationals if there is a positive 5 such that 
\q9—p\ > dq^ 1 holds for all integers p, q . There are uncountably many such numbers; 
nevertheless, they are rather exceptional in the sense that they form a zero-measure 
set. Its complement to the set of all irrationals consists of numbers which we shall call 
Last admissible. 5 A convenient way to characterize these number types is through 
their unique continued-fraction representations: 9 is badly approximable if and 
only if the infinite sequence of integer coefficients in this representation is bounded, 
and Last admissible otherwise. Needless to say, rationals have finitely many nonzero 
coefficients. 

If 9 is irrational, the rhs of (Re fdelta gap) is never zero; the existence of gaps 
requires then that there is a subsequence of local minima which tends sufficiently 
fast to zero. In this way we were able in Ref.13 to prove the following results: 



(a2) 
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(a7) For a badly approximable 9 there is «o > such that for \a\ < «o the 

spectrum has no gaps above the threshold. 
(a8) The number of gaps is infinite for any 9 provided \ot\L > 5 _1 / 2 7r 2 ; recall that 

L : = max(£i, £ 2 ) ■ 

(a9) If 9 is rational or Last admissible, there are infinitely many gaps for any 
a ^ 0. 

The worst irrational in this sense is the golden mean 9 = 1(1 + y/E) which 
has the continued-fraction representation 9 = [1,1,...]. In this case the sufficient 
condition for the existence of infinitely many gaps is necessary at the same time and 
coincides with the critical value of the claim (a7): \cto\L is tt 2 (59)- 1 / 2 = 3.4699... ; 
more about that will be said in Section V below. 



III. 2 5' 3 lattices 

Replacing the S coupling by the boundary conditions (Refdelta' be), one can derive 
the band equation in this case: 13 

cos + cos k£\ cos $ 2 £ 2 + cos k£ 2 (3k 

sin kit + sin k£ 2 2~ = ' ^ 

the same argument as above then shows that spectral bands of the 5' s lattice with 
a coupling constant (3 are determined by the inequalities 

TF T (k)>±^ (6) 

for ±/3 > and k > , and an analogous relation for the negative part. 

The structure of the spectrum is now different; the condition (|6|) allows to make 
the following conclusions: 

(bl) The spectrum equals [0, oo) if and only if /3 = ; otherwise there are infinitely 
many gaps. 

(b2) If (3 > , the lower end of each band coincides with some k 2 or k 2 ^ , where 

n, m are integers. The same is true for (3 < and the upper band ends, with 

the exception of the first one. 
(b3) The lowest band threshold is positive for (3 > and negative if (3 < ; in 

the case —£\ — £ 2 < (3 < the whole first band is negative, and the second 

one starts at zero. 
(b4) The positive bands shrink with increasing \f3\ . 

(b5) Each gap is contained in the intersection of a pair of 5' Kronig-Penney gaps 2 
with the coupling constant (3 and spacings £\ and £ 2 , respectively. 

Instead of the asymptotic behavior (a6) of S- lattices, we have slightly more 
complicated result. If a band high in the spectrum is well separated, its width A r 
is the same as in the 5' Kronig-Penney model, A r = -J- + Oir^ 1 ) . Furthermore, 
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Bild2.1.epsi 

Figure 2: The golden-mean 5 lattice spectrum as a function of the coupling con- 
stant. 



Bild3.epsi 

Figure 3: 5 lattice spectrum as a function of the ratio 9 for a = 20 (above) and 
a = -20 (bottom). 

if 9 is rational and two bands intersect, k n = k m for some n,m , we have a similar 
expression with ij 1 replaced by £i + £ 2 . ^ ma Y happen, however, that k n and 
k m are not identical but close to each other, so that they still produce a single 
band. Then the band width is enhanced; the effect is most profound just before the 
band splits. Using the condition (H), it is straigtforward to estimate the factor of 
enhancement by conspiracy of bands: 13 asymptotically we have the relation 

A + Qir- 1 ) < A r < H (i^+£^) + Oir- 1 ) . (7) 

This brief survey shows, in particular, that despite the generally irregular pat- 
tern, the spectra of 5 and 5' s lattices inherit the main feature of their Kronig-Peney 
analogues, namely that they are dominated asymptotically by bands and gaps, re- 
spectively, at high energies. At the same time, the above listed results leave many 
questions open about the actual form of these spectra and related quantities. In the 
following sections we are going to answer some of them. 

IV 5 lattice spectra 

The form (£|) of the gap condition makes it possible to find the 5 lattice spectra 
numerically in dependence on parameters of the model; the results are shown on 
Figures 2 and 3. The effect of competing periods is obvious: while a square lattice 
has a familiar Kronig-Penney spectrum shape 2 with the halved coupling constant 
as expected from the gap condition, in the general case the gap pattern is irregular. 

The dependence of spectral bands on the ratio 9 shows clearly that the gaps 
are contained in view of (a5) in the intersections of KP-gaps. The latter can be 
numbered by an integer n such that nnij 1 is the lower (upper) gap end for ±a > , 
respectively (with n = referring to the region below the bottom of the spectrum 2 ). 
This allows us to label the gaps of a S lattice naturally by a pair (ni, n 2 ) with the 
integers rij tagging the KP-gaps. Notice also that due to (a2) the actual gaps 
occupy the lower part of the allowed region if a > , and the upper one for a < . 
Asymptotically each intersection is of a parallelepipped form and the actual gap 
part represents a diagonal cut of it, which conforms with (a6). 
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Bild4. l.epsi Figure 4a: The folded 5 lattice spectrum as a function of ln^: The 
regions with at least one gap; the lowest one, (1,1), is deleted. 



Bild4.2.epsi Figure 4b: The folded 5 lattice spectrum as a function of In 9: An 
inset of the previous picture. 

Although the 9 plots look rather regularly, they contain a hidden irregular 
pattern which is revealed if the spectrum is properly folded. A way to do that 
is suggested by a relation between the present model and the theory of multidi- 
mensional Jacobi matrices 21 by which A = 2 sm(k(£i + £ 2 )) is the natural "folded" 
energy variable, the above introduced gap labelling has then a direct relation to that 
of Ref.3. The folded spectrum is illustrated on Figure 4; we plot A against In 9 to 
show the symmetry with respect to the exchange of the basic cell sides. All the gaps 
close, of course, and the multiplicity is infinite, but it changes in a fractal way 22 
as it is seen when we plot the regions where the original spectrum has at least one 
gap. For the sake of illustrativeness, we delete here the "regular" (1, 1) gap which 
produces a concave strip covering a part of the upper half of the picture. Similar 
"non-fractal" contributions come from the gaps [n%, 1) and (l,n2) with rij > 2 
which enter the picture at sufficiently large | ln#| . Each gap of the fractal structure 
"contains" an infinite series of embedded gaps as can be seen if we plot just the gap 
edges; for instance, the central element (2, 2) of the picture is, in fact, a sequence 
(2,2)D(3,3)D(4,4)D.... 

V Critical coupling constants 

Due to (a7), the spectrum may contain no gaps if the ratio 9 is badly approximable 
and the coupling is weak enough. To get a better understanding of how the number- 
theoretic properties of the parameter imply this result, let us discuss in more detail 
the "worst" case of the golden-mean lattice, 9 = |(1 + a/5) . 

Suppose that the coupling constant grows from the zero value. The condition 
(|4j) shows that a new gap opens whenever the hyperbolic graph of the function on 
the Ihs crosses a local minimum of the function F + (-) . These critical values of the 
coupling together with the corresponding momenta are shown on Figure 5. One can 
make from here several conclusions: 



Bild4.3.epsi Figure 4c: The folded 5 lattice spectrum as a function of ln#: The 
plot of gap edges in the same scale as (a). 
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Bild4.4.epsi 

Figure 4d: The folded 5 lattice spectrum as a function of In 9 : Gap edges in a 
wider scale, with several lowest of the (ni, 1) and (l,n 2 ) gaps shown. 

Bild5.epsi 

Figure 5: Critical coupling constants for the golden-mean 5 lattice. 

(i) Gaps open in series. The asymptotic behavior is clearly visible and the critical 

points approach the asymptotical value from above; hence a new infinite series 
of gaps opens "from above" once the coupling constant crosses the next critical 
value. 

(ii) The critical points are extremely sparse; we see that the patterns are practically 

equidistant in the logarithmic scale. This is a nice illustration of the fact that 
9 is badly approximable: to obtain the next close approximation to it we have 
to use integers which are roughly twice as large as their predecessors. 

(hi) The asymptotic critical values are quantized; they appear at multiples of the 
two basic values Z/g ±1//2 ) where I := \/t\l2 ■ However, not every multiple 
yields a critical value: in both series the same sequence of integers repeats, 
namely 

1,4,5,9,11,16,19,20,25,... (8) 

Let us attempt to explain this quantization rule. Recall that the golden mean is 
approximated by ratios of successive Fibonacci numbers, 17 9 = lim^oo Un+iu^ 1 , 
where the u n satisfy the recursive relation u n+ i = u n + u n ^i and assume the values 

1, 2, 3, 5, 8, 13, 21, . . . for n = 1,2, In the continued-fraction representation the 

n-th approximant equals 9 n = [1, . . . 1, 0, . . .] with n + 1 nonzero coefficients. 

Let us "spoil" the approximation by choosing instead the sequence of 0% '■ = 
[1, . . . , 1, N, 0, . . .] with a fixed positive integer N preceded by n ones. This can 
be alternatively written as gn+ig" 1 , where q n := Nu n + u n -\ , so that 



2 | Q 'I" ■ 1 1 



Vn J 



N 9 n - {-9)- n n - 1 -(-0)~ n + 1 



-1)' 



1 + N - N 2 



N 
J 



V5 

as n — > oo . The asymptotic threshold values thus are 

r 2 



n ' #±1/2 



N 2 -N-l 
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which yields the integer multiples 1,5,11,19,29,41,... explaining a part of the 
values (||) but not all of them. The missing numbers can be obtained if we consider 
a more general approximation to the golden mean, 

qA.b Au n+2 + Bu n+ i 
n ' Au n+1 +Bu n 

for some positive A, B . The numerical factor in the above limit is then replaced by 
\A 2 — B 2 — AB\ , so the complete sequence ([8]) is obtained when A,B run through 
positive integers. What is important is that A, B here need not be relatively prime. 
This can be seen directly: it is clear from (|J) that replacing approximation sequences 
{in}, {Pn} by {Mq n }, {Mp n } means that the corresponding local minima of F + {k) 
are enhanced by a factor which tends to M as n — > oo . These are compared with 
at the minima which is divided by M , so the asymptotic critical value of the 
modified sequence is multiplied by M 2 . In the logarithmic scale the x coordinates 
referring to the modified sequence are just shifted by logM ; this is also seen on 
Figure 5. 

In this way, we are able to explain how the successive infinite gap series open in 
the spectrum as the coupling constant a increases. 



VI Fermi surfaces 

The above results do not exhaust all possible irregularities of the model in question. 
The fact that a given energy value is contained in a spectral band tells us nothing 
about the directions in which the electron may propagate. To this end we have to 
find what are the allowed values of the quasimomentum, i.e., to determine the Fermi 
surface. As mentioned above, the Brillouin zone is the rectangle 
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(10) 



and the Fermi surface is a one-dimensional submanifold in it. Since it has a four-fold 
symmetry due to the invarince of the condition ([$]) with respect to the interchanges 
dj —dj , it is sufficient to consider a quadrant of B . The change of variables 
Vj := cos$ j£j maps it on the square [—1, 1] x [—1, 1] ; the band condition becomes 
at that 

+ ,„ = — + cot kit + cot k£ 2 



sin k£\ sin k£ 2 Ik 

describing a line segment which crosses the square provided k corresponds to a 
point in a spectral band. Passing back to the quasimomentum coordinates we see 
that the Fermi surface can have one of several standard shapes. 

On the other hand, the way in which the surface chang function of the 

Fermi energy depends substantially on the ratio 9 . This is illustrated on Figure 6 
where we plot the two-dimensional manifolds spanned by the Fermi surface as the 
energy increases; to see the change over a larger range of energies we choose for the 
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Bild6.1.epsi Figure 6a: Fermi surfaces for different energies for the second band of 
a square 5 lattice. 



Bild6.3.epsi Figure 6b: Fermi surfaces for different energies, the "side view" for the 
same lattice and the logarithmic scale of energy, 1 < k < 21, . 

latter a logarithmic scale. We see that for a square lattice the surface "evolution" 
produces a regular pattern of switching "caps" reminiscent of the corresponding 
quantity for the one-dimensional Kronig-Penney model 2 , while in the golden-mean 
case the dependence is highly irregular. 

Another lesson from the dependence of Fermi surfaces on energy concerns the 
character of the spectrum. To avoid spectral singularities, the solution k = k(9\, 82) 
of the equation @ has to be smooth with local extrema at discrete points only. The 
Ihs of (Refdelta band condition) which we denote as D(k, 9%, 62) is certainly smooth 
in all variables away of the points k£j = wk, and 

l? (MlA) = ^FMl + siJk + 2^ ^ 

holds there, with a possible exception of a discrete set of points for a < . Hence 
k(0i, 82) is smooth by the implicit-function theorem and 

dk (a ft \ - p sin9 ^ ( dD (h ft 

which may be simultaneously zero only if 8j£j = ±7r , j = 1,2. This means that 
stationary points of k(8i,8 2 ) are only at the center of the Brillouin zone or at its 
cornerpoints; we conclude that 

(alO) the 5 lattice spectrum in the bands is absolutely continuous. 

This is not just a mathematical statement; it tells us that the transport properties 
of electrons in the allowed energy windows are "normal" . 

VII External fields 

Adding an external field, even a simple one like a linear potential, makes the prob- 
lem much more difficult, and we are not able to present more than several heuristic 




Bild6.2.epsi Figure 6c: Fermi surfaces for different energies for the second band of 
a golden-mean S lattice 



11 



Bild6.4.epsi 

Figure 6d: Fermi surfaces for different energies, the "side view" for the same lattice 
and the logarithmic scale of energy, 1 < k < 21, . 



observations. Recall that the problem is far to be understood even in the one- 
dimensional case. It has been demonstrated recently that an array of 5' interactions 
has no absolutely continuous spectrum when an electric field is applied; 10,23 this 
result extends to a wide class of background potentials. 24 Moreover, one may con- 
jecture that the spectrum of 8' Wannier-Stark ladders depends substantially on the 
slope of the potential being pure point and nowhere dense if the latter is rational, 
and covering the whole real axis otherwise. 23 

The case of a periodic 5 array is more difficult and only partial results are 
known, 25-27 though it is generally believed that a phase transition occurs in this 
case with the spectrum being pure point at weak fields and continuous if the field is 
strong. 26 The critical field value can be conjectured to be F crit « a 2 /4a , where a is 
the array spacing and F the field intensity, using the argument of Ref.26, which is 
based on estimating the tunneling probability through the family of tilted spectral 
gaps. 28 The result is intuitively appealing because it compares two value having an 
established meaning: the energy step Fa with the bound-state energy a 2 / 4 of a 
single 5 well. However, the argument neglects the fact that in some tilted gaps the 
other (exponentially growing) solution for the classically forbidden region may play 
role, and therefore it is not apriori clear whether this is the correct answer. 

Applying the same tilted-band picture to the two-dimensional situation, one 
may conjecture that in a S' s lattice whose spectrum is dominated by the gaps at 
high energies, an external electric field is likely to produce a localization in the field 
direction; in other words, electrons would be able to move at most in the direction 
perpendicular to the field. 

In the S lattice case the situation is more complicated. Even if a Berezhkovski- 
Ovchinnikov-type argument 26 could be justified, it will yield now only a lower bound 
to the power with which the solution to the Schrodinger equation decays, because a 
gap-width counting disregards the fact that at some energies belonging to a band a 
propagation is possible in certain directions only. With this reservation in mind, let 
us compute the expression which is a two-dimensional analogue of the BO-power: 
it equals (n/16F)T(9,a) with 



where 5 m = 5 m (9,a) are the corresponding gap widths. 

Due to the results of Section IV, namely the property (a4), the sum in ([11]) is 
monotonously increasing as a function of \a \ , the growth being roughly quadratic. 
The same behavior is expected for T(8, a) though for a badly approximable 9 the 




(11) 
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Bild7. l.epsi 

Figure 7: The quantity T(9, a) as a function of the coupling constant for 9 = 
1, |, HI , and the golden mean. The arrows mark distinguishable gap openings. 



Bild8.3.epsi Figure 8a: The quantity T(9, a) as a function of 9 . Plot for 
a = 3.4699 over rational values of 9 . 

existence of the thresholds discussed in Section 5 may lead to local deformations of 
this dependence around the threshold values of a ; this is illustrated on Figure 7. 

On the other hand, let us fix the coupling constant at two values, a = 3.4699... 
for which we know that T{9, a) has a zero, and a = 20 , and plot the 9 depen- 
dence (Figure 8). Simple harmonic values of the ratio produce well pronounced 
peaks; hence the corresponding lattices are expected to exhibit better localization 
properties when an electric field is applied. On the other hand, T(p/q,a) may 
depend substantially on the integers p, q involved. To illustrate the difference we 
compute first this quantity for p, q = 1, 2, 3, . . . with q < 150 , and after that we 
plot T(6>, a) with random identically distributed 9 values; the latter are almost 
surely "typical", i.e., irrational. 

VIII S f s lattice spectra 

Let us finally check how the spectra of S' s lattices look like. They can be found with 
the help of the condition (^|); in the same way as above, we can also check that 

(b6) the S' s lattice spectrum in the bands is absolutely continuous. 

Figure 9 shows the dependence of spectral bands on (3 , again for the golden-mean 
case, 9 = |(l+v2) ; we see that the pattern has globally the same behavior as the 5' 
Kronig-Penney spectrum. 2 On the other hand the 9 plot (Figure 10) illustrates the 
properties (b2) and (b5); the latter can be again used to label the gaps. Moreover, 
we see clearly the enhancement due to conspiracy of bands - cf. (0) . 

As in the S situation, the spectrum may be folded into 2 sin(/c(£i + £ 2 )) ; this 
is shown on Figure 11. We have deleted here again the lowest "regular" band 
which can be done unambigously if \ln9\ is small enough. In this case the folded 
spectrum leaves, in general, some gaps open, but it may happen at rational values 
of 9 only. Indeed, it is not difficult to check that for an irrational 9 spectrum 



Bild8.4.epsi Figure 8b: The quantity T(9, a) as a function of 9 . Plot for 
a = 3.4699 with 9 being random identically distributed numbers. 
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Bild8.1.epsi Figure 8c: The quantity T(9,a) as a function of 9 . Plot for a = 20 
over rational values of 9 . 



Bild8.2.epsi 

Figure 8d: The quantity T(9, a) as a function of 9 . Plot for a = 20 with 9 being 
random identically distributed numbers. 



Bild9.2.epsi 

Figure 9: The coupling constant plot for the S' s lattice spectrum in the golden mean 



case. 



Bildl0.2.epsi 

Figure 10: The S' s lattice spectrum as a function of the 9 for (3 = 5. 



Bildll.epsi 

Figure 11: The folded 5' s lattice spectrum as a function of the 9 for a = 20 . The 
value of momentum is restricted to k < 100 . 



Bildl2.1.epsi Figure 12a: Folded-spectrum gaps for rational approximations to the 
golden mean: (3 = 20, . 



Bildl2.2.epsi 

Figure 12b: Folded-spectrum gaps for rational approximations to the golden mean: 
[3 = 5. 



14 



the folded spectrum has a full measure. It is sufficient to notice that all the points 
2? 9 ±1/2 with both signs are contained in a band, either as its lower edge (upper 
edge for (3 < ) or an internal point. The corresponding values of A are then 
2 sin(7m(l + # ±1 )) ; they cover the interval [—2,2] due to the basic result of the 
ergodic theory. 29 On the other hand, for rational 9 some gaps remain open; this is 
illustrated on Figure 12. 

IX Conclusions 

We have demonstrated irregular properties of rectangular 5 and 5' s lattices with 
incommensurate sides, as they are manifested in the form of their spectra, critical 
coupling-constant values, dependence of the Fermi surface on the Fermi energy, and 
possibly in localization properties under influence of external fields. These features 
are expected to play a dominating role in behavior of finite but large enough graph 
superlattices. It is also natural to conjectures that similar geometrically induced 
spectral properties will be seen in superlattices in which the basic cell has another 
form. 

One lesson drawn from the present considerations concerns the validity of the 
so-called Bethe-Sommerfeld conjecture by which the number of gaps in periodic 
systems of dimension two or higher is finite. This was proven by SkriganovQ for 
Schrodinger operators with nice potentials, and it is known to be valid also for some 
systems with S interactions. 2 We have seen that the conjecture may not be valid if 
the periodic structure is singular enough: the golden-mean 5 lattice has, depending 
on the value of the coupling constant a , either infinitely many gaps or none at all 
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